Collective stimulated Brillouin backscatter 
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We develop the statistical theory of the stimulated Brillouin backscatter (BSBS) instability of 
a spatially and temporally partially incoherent laser beam for laser fusion relevant plasma. We 
find a new regime of BSBS which has a much larger threshold than the classical threshold of a 
coherent beam in long-scale-length laser fusion plasma. Instability is collective because it does not 
depend on the dynamics of isolated speckles of laser intensity, but rather depends on averaged 
beam intensity. We identify convective and absolute instability regimes. Well above the incoherent 
threshold the coherent instability growth rate is recovered. The threshold of convective instability 
is inside the typical parameter region of National Ignition Facility (NIF) designs although current 
NIF bandwidth is not large enough to insure dominance of collective instability and suggests lower 
instability threshold due to speckle contribution. In contrast, we estimate that the bandwidth of 
KrF-laser-based fusion systems would be large enough. 

PACS numbers: 52.38.-r 52.38.Bv 



Inertial confinement fusion (ICF) experiments require 
propagation of intense laser light through underdense 
plasma subject to laser-plasma instabilities which can be 
deleterious for achievement of thermonuclear target igni- 
tion because they can cause the loss of target symmetry 
and hot electron production ^ . Among laser-plasma in- 
stabilities, the backward stimulated Brillouin backscat- 
ter (BSBS) has been considered for a long time as seri- 
ous danger because the damping threshold of BSBS of 
coherent laser beams is typically several order of mag- 
nitude lower compared with the required laser intensity 
^ lO^^W/cm^ for ICF. Recent experiments for a first 
time achieved conditions of fusion plasma and indeed 
demonstrated that large levels of BSBS (up to tens per- 
cent of reflectivity) are possible Q- 

Theory of laser-plasma interaction (LPI) instabilities 
is well developed for coherent laser beam [3|. However, 
ICF laser beams are not coherent because temporal and 
spatial beam smoothing techniques are currently used to 
produce laser beams with short correlation time, Tc, and 
lengths to suppress laser-plasma interactions. The laser 
intensity forms a speckle field - a random in space dis- 
tribution of intensity with transverse correlation length 
Ic — 2F/ko and longitudinal correlation length (speckle 
length) Lspeckie — 7F^Ao, where F is the optic //# and 
Ao = 27r/fco is the wavelength (see e.g. 0, Q)- Beam 
smoothing is a part of most constructed and suggested 
ICF facilities. However, instability theory of smoothed 
laser beam interaction with plasma is not well developed. 
There are intense experimental and simulation ongoing 
efforts 0] to determine BSBS threshold for smoothed 
beams which appears to be in some cases quite low so 
that it is now under discussion that laser intensity at the 
National Ignition Facility (NIF) should lowered by a fac- 
tor of few compared with original NIF designs [l| with 
intensities ~ 2 x lO^^W/cm^. 

Here we develop a theory of collective BSBS instabil- 



FIG. 1: Regimes of BSBS of partially incoherent laser beam. 
Tc is the laser correlation time and jbsbs = ijJi is the BSBS 
temporal grow rate, both in arbitrary units. For large Tc in- 
stability is dominated by speckles of laser field (on the right). 
CBSBS regime is shown in the middle. Very small Tc corre- 
sponds to RPA regime (on the left). Solid lines sketch tran- 
sitions between all 3 regimes. 



ity (CBSBS), which is a new BSBS regime, for prop- 
agation of laser beam with finite Tc in homogeneous 
plasma. CBSBS has threshold comparable with NIF in- 
tensities. CBSBS requires T^ small enough to suppress 
contribution from speckles. If we additionally assume 
that Tc 3> L speckle /c then CBSBS threshold does not 
depend on Tc. Such Tc is accessible to KrF lasers [(i], 
Tc ~ 0.7ps, but not for NIF glass lasers with beam 
smoothing up to 3 A at Ilu, implying Tc ~ 4ps at 3uj. 
This is consistent with the numerical simulations which 
show that BSBS threshold in NIF emulation experiments 
is dominated by speckles We show below that 

speckle-dominated threshold is lower by a factor 7 than 
CBSBS threshold. Since plasma inhomogeneity can only 
increase instabihty threshold The CBSBS threshold 
is a lower bound. Fig. [T] depicts CBSBS between large 
Tc speckle regime 0] and random phase approximation 
(RPA) regime. 

Assume that laser beam propagates in plasma with 
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frequency ujq along z with the electric field £ given by 
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where E{r,z,t) is the envelope of laser beam and 
B(r,z,t) is the envelope of backscattered wave, r = 
{x,y), and c.c. means complex conjugated terms. Fre- 
quency shift Alu = —2kQCs is determined by coupling of 
E and B through ion-acoustic wave with phase speed Cs 
and wavevector 2ko with plasma density fluctuation Sug 
given by ^ = itTe^*'^"^+''^"*-|-c.c., where a{r, z,t) is the 
slow envelope and rie is the average electron density, as- 
sumed to be small compared to critical density, ric. The 
coupling of E and B to plasma density fluctuations gives, 
ignoring light wave damping. 
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V — {dx,dy), and a is described by the acoustic wave 
equation coupled to the pondermotive force oc £^ which 
results in the envelope equation 

[i{c;'dt + 2v,ako + dz) - (4fco)~' V2]a* = ~2koE*B, (4) 

where we neglected terms oc \E\'^, \B\'^ in r.h.s. which 
are responsible for self-focusing effects, vl is the Lan- 
dau damping of ion-acoustic wave and = VL/2koCs is 
the scaled acoustic damping coefficient. E and B are in 
thermal units (see e.g. |12|). 

Assume that laser beam was made partially incoher- 
ent through induced spacial incoherence beam smooth- 
ing [l^ which defines stochastic boundary conditions at 
z = for the spacial Fourier transform (over r) compo- 
nents -E(k), of laser beam amplitude 12]: 

E{k,z = 0,t) = \Ei,\exp[i<j)i,it)], 
(exp i [0k (t) - (p^, {t')]) = <5kk' eM-\t-t'\/Tc), 



which means that we neglect off-diagonal terms 
-^kj'i j 7^ /• Since speckles of laser field arise from 
interference of different Fourier modes, j ^ j', wc as- 
sociate the off-diagonal terms with speckle contribution 
to BSBS (independent hot spot model [1, ) • Speckle 
contribution can can be neglected if [l6l | 



(7) 



where tsat is the characteristic time scale at which BSBS 
convective gain saturates. 

We use the linear part of the theory of Ref. 
to estimate T^at for speckle contribution to backscat- 
ter as tsat = {L speckle / c)[2 + (7o/j^l)^], where 7^ = 
kQCCsIspeckiene/2nc and we choose the typical intensity 
of light in speckle Ispeckie = 3/, where / is the spatial 
average of laser intensity \E\'^ [1, [l^, [3]. In such a case 
Tc/tsat — TekgCsVia / (4,1) , where here and below / desig- 
nates the scaled dimensionless laser intensity defined as 
i = ^^I. For typical NIF parameters / ~ 1 [H, [13, 
Aq — 351nm and Cs — 6 x 10^ cm we obtain from 
([7]) the estimate ^ OA/via [ps] which is not satis- 
fied for low plasma ionization number Z plasma in NIF 
which typically has I'ia ^ 0.1. However, CBSBS can 
still be relevant for NIF in gold plasma near hohlraum 
^ with Via ~ 0.01. Similar estimate for KrF lasers 
(Ao = 248nm, F = 8) gives Te < O-S/via [ps] which is 
easier to satisfy because of smaller Tc and suggests that 
KrF lasers are better suited for applicability of CBSBS. 

If we look for solution of Eqs. ([3]) and ^ in exponen- 
tial form Bj,(T* oc e'lf^^+k r-ajt)^ arrive at the follow- 
ing dispersion relation in dimensionless units 



iuj + fi + IK — [i/A)k'^ 



N 

J = l c 
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where /i = 2viak'^/k'^, I /km is the transverse unit of 



£'k| = const, k < km; -Ek = 0, k > km, (5) length, ko/k^ is the unit in z direction, ko/k'^Cs is the 



chosen to the idealized "top hat" model of NIF optics 
14| . Here km — kQ/{2F) and the average intensity, (/) = 
(|£'|^) = / determines the constant. 

In linear approximation, assuming ^ \E\ so that 
only laser beam is BSBS unstable, we can neglect right 
hand side (r.h.s.) of Eq. The resulting linear equa- 
tion with top hat boundary condition ([S]) has the ex- 
act solution as decomposition of E into Fourier series, 
E{r, z, t) = Y,j Ek- with E^j oc exp [i{(j)kj [t -z/c) + kj ■ 
r — k|z/2fco)] . Eq. ([3]) is linear in B and E which implies 
that B can be also decomposed into B = B\s_. . We 
approximate r.h.s. of @ as E*B ~ E^,B\^^ so that 

-1^ 



= -2koY^E^ 



Bu 



(6) 



time unit and / = J2j \Ekj 

The dispersion relation ([5]) is correct provided the tem- 
poral growth rate Ui — Im{uj) is small compare to inverse 
time of light propagation along speckle, cjj <C c/ LspeUe, 
and if during time Tc light travels much further than a 
speckle length, L speckle cTc- That second condition 
ensures that term oc (f)'^\t — z /c) ~ 1 /Tc could be ne- 
glected in Eq. ([3]) allowing the time dependence of E 
in Eqs. ([3|) and ([6|) to be ignored and in such case den- 
sity fiuctuation a evolves without fiuctuations. E.g. for 
typical NIF parameters, TckoCs/2F 1 we obtain that 
2c/7csF > 1 which is well satisfied for NIF optics [L]. 

In the continuous limit N ~* 00, sum in ([5]) is replaced 
by integral which gives for most unstable mode k = 0: 

A(uj,k) ^ ~iuj + ri + iK + i^iln = 0. (9) 

4 -K - UJ^ 
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Eq. ([9]) has branch cut in complex k plane determined 
by points ki = 1 — lo^ and K2 = —f^^- Standard anal- 
ysis of convective vs. absolute instabilities (see e.g. fisj l 
should be modified to include that branch cut. In dis- 
crete case with N ^ 1 instead of branch cut the dis- 
crete dispersion relation ^ has solutions located near 
the line (/ci,K2). These solutions are highly localized 
around some so they cannot be approximated by ([9|) 
but they are stable for iV ^ 1. Generally there are two 
solutions of ([9]), however for Im{uj) oo one solution 
is absorbed into branch cut. Second solution is stable. 
Above the convective CBSBS threshold, 

Iconvthresh = 4/7r, (10) 

the first solution crosses real n axis from below as 
Im{u!) ^ so it describes instability of backscattered 
wave with Im{K) > 0. 

However, above the absolute CBSBS threshold, which 
can be approximated from solution of Eq. ^ as 

Lbsthresh ^ (1/2) (a*"^ + a* + V/^^ - 2) , fi>4, (11) 

the contour Im{uj) — Const cannot be moved down to 
real lo axis because of pinching of two solutions of ([9]) 
which defines growth rate of absolute instability. We con- 
clude that classical analysis of instabilities still holds for 
incoherent beam if we additionally allow the absorption 
of one solution branch into branch cut. This effect re- 
sults from incoherence of pump beam which has infinitely 
many transverse Fourier modes in approximation of Eq. 
([9|) and there is no counterpart of that effect for coherent 
beam. 

For /i ^ 1 the absolute threshold pT|) reduces to the 
coherent absolute BSBS instability threshold 

-^absthreshcoherent — (12) 

For NIF parameters, ~ 5keV, F = 8, Ue/ric = 
0.1, Ao — 351nm with moderate acoustic damping, ~ 
0.1, we obtain in dimensional units Iconvthresh — 2 x 
lO^^W/cm^ and labsthresh ^ 9 X IQi^W/cm^. For high Z 

asma (e.g. gold plasma near the wall of NIF hohlraum 

, v^a ^ 0.01) we obtain Iconvthresh ^ 2 X lO^^'W/cm^ 
and labsthresh ^ 9x lO'^'^W/cm^. Typical intensity of NIF 
laser shots is between lO^^W/cm^ and 2 x lO^^W/cm^ so 
we conclude that in different parts of NIF plasma both 
convective and absolute instabilities are possible. Fig. [2] 
compares instability gain rate of coherent and incoherent 
beams for fi = 51.2. 

In contrast with Eq. pO|l . the convective instability 
threshold in coherent case is because we neglect damp- 
ing of B in Eq. ([3]). Retaining coUisional light damping 
gives finite threshold Iconvcoherent = IdF'^i'B/koc < 1, 
where i^b = [3| is the coUisional damping of 

backscattered wave B and Vci is the electron-ion colli- 
sion frequency. That threshold is several orders of mag- 
nitude smaller compared with (jlOp and is neglected here. 
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FIG. 2: Convective instability gain rate Im{K) as a function 
of dimensionless laser beam intensity (lower grid) for inco- 
herent laser beam (solid curve) and coherent beam (dashed 
curve) for fi — 51.2. Here n = SF'^Uia- Upper grid corre- 
sponds to laser intensity in dimensional units for typical NIF 
parameters Te ~ 5keV, F = 8, Ue/ric — 0.1, fia = 0.1, ujo ~ 
3.6 X lO^^sec"^. Note that in the range between Via = 0.05 
and Via = 0.3 there is no qualitative change in the behavior 
of both curves except change of scale at upper axis. 
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FIG. 3: Angular width 2FAk/ko of convective instability gain 
rate Im^n) vs. laser intensity scaled to (|10|l for /i = 153.6 
There is no qualitative change in that curve for smaller fj, 
down to /i = 5.12. For NIF optics, F = 8, value p — 153.6 
corresponds to Via = 0.3 (e.g. He-H plasma) and = 5.12 
corresponds to Via ~ 0.01 (e.g. gold plasma). 

Qualitatively incoherence of laser beam can be consid- 
ered as effective damping of B with effective damping 

rate ineffective — \QF^ ' 

Depending on laser incoherence we have a hierarchy of 
thresholds: (a)Spatially incoherent laser beam with large 

Tc has threshold, Ithreshold = Ithresholdsveckle = 4/77r 

which is dominated by intense speckles |7|]. (b) Spatial 
and temporary incoherent beam with Tc satisfying ([7]) is 
given by ()10|) which factor 7 times higher compared with 
speckle threshold and does not depends on Tc- It indicates 
practical limit of how threshold of BSBS instability can 
be increased by decreasing Tc. (c)For much smaller Tc, 
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such that it is smaller than both inverse acoustic damping 
Tc <C and inverse temporal growth rate Tj. <C l/w^, 
the classical RPA regime is recovered which has ignorable 
diffraction ([i,[ll[ll|). This limit (e.g. for Aq = 351nm 
and Via = 0.15 it requires Tc <C 0.3ps) is not practical for 
ICF as Tc is too small. Cases (a)-(c) are shown in Fig. 

HI 

Current NIP 3A beam smoothing design is between 
regimes (a) and (b). KrF laser with Tc ~ 0.7ps would 
be in regime (b). Thus generally we expect that KrF- 
laser-based ICF allows access to CBSBS regime although 
CBSBS threshold for NIF can be possibly initiated by 
self-induced temporal incoherence (see e.g. [3]). An- 
other possibility for self-induced temporal incoherence 
is through collective forward stimulated Brillouin scat- 
ter (CFSBS) instability [il,^. Above CFSBS thresh- 
old correlation length decreases with beam propagation 
length and may decrease Tc. For low Z plasma threshold 
for CFSBS is close to (fTO|) 12]. As Z increases (which can 
be achieved by adding high Z dopant), CFSBS threshold 
decreases below ([T0| and might result in decrease of Tc- 

To distinguish contribution to BSBS from speckles 
(regime (a)) and CBSBS (regime (b)) we propose to look 
at angular divergence /\6 — Ak/ko of BSBS. In gen- 
eral one expects gain narrowing of the scattered light: 
the modes close to the most unstable mode, with gain 
rate {Ki)max, dominate. Here Ki = Im^n). Fig. [3] 
shows 2FA9 from CBSBS as a function of laser in- 
tensity above CBSBS threshold at propagation distance 
L — 10/{Ki)„iax- L is chosen from the physical condition 
that there is sufficient convective CBSBS gain, to amplify 
the energy of thermal acoustic fluctuations at wavenum- 
ber 2fco to have reflectivity ^ 1, and for fusion plasma 
this is typically exp(G) = exp(20) (see e.g. [3]), where 
G = 2KiL is the power convective gain exponent. Then 
Ad is conventionally defined by half width at half max- 
imum: exp[G(A0)] = O.5exp[G(0)]|e=o. Important fea- 
ture of CBSBS seen in Fig. El is that A6' ^ at threshold 
with Ki{k) ~ K,;(0)(1 — ak^) and a ~ fil /{^.I — 1) near 
threshold. Fig. [3] should be compared with AO from 
speckle-dominated backscatter. Previous work QQ\ sug- 
gested that speckles can also cause A 9 below top- hat 
width, 1/2F, for very intense speckle backscatter. We 
estimate based on Refs. 0, i]) that for nominal ICF 
plasma (^ 10^ speckle volumes), most intense speckle is 
~ 15/ which gives Gmtense — 15(Go) — 100 near CBSBS 
threshold, where (Gq) = 2tiiL speckle — 7 is the the gain 
over speckle with the average intensity /. We performed 
direct simulations of backscatter from Gintense — 100 
speckle and found that A^ ~ 1/2F which means that 
asymptotic [i^l is still not applicable. In other words, fi- 
nite size plasma effects dominates over asymptotic theory 
of inflnite plasma. We conclude that regime (a) can be 
easily distinguished from CBSBS regime (b): near CB- 
SBS threshold with condition ([7]) satisfied one should see 



backscattered light spectrum with essential peak whose 
width is given by Fig. [3] and wide weak background de- 
termined by speckles. 

In summary, we found a novel coherence time regime 
in which Tc is too large for applicability of well-known 
statistical theories (RPA) but rather an intermediate 
regime, Tc is small enough to suppress speckle BSBS. Un- 
like coherent beam CBSBS has threshold typically much 
larger than that determined by damping while for laser 
intensity many times above convective instability thresh- 
old for incoherent beam, the coherent theory is recovered. 
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